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Abstract 

We report studies of semileptonic decays, B — ► X c £u, based on a sample of 88 million BB events 



recorded with the BABAR detector. We have measured four moments of the electron energy 
distribution and four moments of the hadronic mass distribution, each as a function of the minimum 
electron energy. From these moments we determine the inclusive branching fraction, the CKM 
matrix element \V c b\-, and other heavy quark parameters, using Heavy Quark Expansions (HQE) 
to order l/rrc| in the kinetic mass scheme. In addition, we have studied a large sample of exclusive 
B° — ► D*~l + v decays. This sample is used to extract the vector and axial form factors, the 
normalization and slope of the HQET form factor to determine \Va\. 
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1 Introduction 



The CKM matrix element V c b, the dominant coupling of the b quark to the charged weak current, 
is one of the fundamental parameters of the Standard Model. It determines the rate for b — ► civ 
decays which at the parton level can be calculated accurately. This rate is proportional to \V c b\ 2 
and depends also on the quark masses, nib and m c . 

2 Inclusive Measurements 

To relate measurements of the inclusive semileptonic l?-meson decay rate to \V c b\, the parton- 
level calculations have to be corrected for effects of strong interactions. Heavy-Quark Expansions 
(HQEs) [1] have become a useful tool for calculating these perturbative and non-perturbative QCD 
corrections [2] and for estimating their uncertainties. We have chosen the kinetic-mass scheme [3, 4] 
for these expansions in l/m^ and a s (m,b) (the strong coupling constant). To order 0(l/m\) there 
are six parameters: the running kinetic masses of the b and c quarks, mb(p) and m c (p), and 
four non-perturbative parameters. The leading non-perturbative effects arise at 0(l/mf) and are 
parameterized by p 2 (p) and Pq(p), the expectation values of the kinetic and chromomagnetic 
dimension-five operators. At 0(l/mf), two additional parameters enter, p%(p) and p\g(p), the ex- 
pectation values of the Darwin (D) and spin-orbit (LS) dimension-six operators. These parameters 
depend on the scale p that separates short-distance from long-distance QCD effects. 

We determine these HQE parameters from a fit to the moments of the hadronic-mass and 
electron-energy distributions in semileptonic B decays to charm particles, B — > X c £v, averaged 
over charged and neutral B mesons. We have measured these moments as functions of E cut , a 
lower limit on the lepton energy (for energy moments we only use electrons, for mass moments we 
also use muons). The moments are corrected for detector effects and QED radiation. The charmless 
contributions are subtracted, using the branching fraction B U £ U =(0.22 ± 0.05)% [5]. 

The hadronic-mass distribution is measured in events tagged by the fully reconstructed hadronic 
decay of the second B meson [6]. The hadronic-mass moments are defined as M^(E cu t) = 
(nix) E e > E cut with n = 1,2,3,4. The electron-energy distribution is measured in events tagged 
by a high-momentum electron from the second B meson [7]. We define the first energy mo- 
ment as Mf(E cut ) = {Ee)Et>E cu t an d the second and third moments as M^(E cut ) = ((Eg — 
Mi(E cut )) n ) E £ >E cut w ith n = 2,3. In addition, we use the partial branching fraction MQ(E cut ) = 
fg max (dBtfv/dEg) dEg. All measured moments and the results of the least-square fit are shown 
in Fig. 1. The fit describes the data well with \ 2 = 15.0 for 20 degrees of freedom. For \V c b\, 
the semileptonic branching fraction, and the heavy-quark masses (at the scale of p = 1 GeV ) we 
obtain [8] 

\V cb \ = (41.4 ± 0A exp ± 0.75 HQE ) nr 3 
B cev = (10.61 ± 0.16 e:cp ± OMhqe) % 

m b = (4.61 ± 0.05 exp ± 0.04^ QE ) GeV 

m c = (1.18 ± 0.07 exp ± 0.06hqe) GeV , 

and m,b — m c = (3.436 ± 0.025 exp ± 0.018hqe ± 0.010 as ) GeV . Beyond the statistical, systematic 
and HQE uncertainties that are included in the fit, the limited knowledge of the expression for 
the decay rate, including various perturbative corrections and higher-order non-perturbative cor- 
rections, introduces an additional error on \V c b\, assessed to be 1.5% [4] and included in the stated 
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HQE error. The choice of the scale (i is estimated to have a very small impact on |V^| and the 
branching fraction [4]. For the non-perturbative parameters in the kinetic-mass scheme up to order 
(1/mf) we obtain 
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Pg 
Pi 
Pls 



0.45 ± 0.04 exp ± O.OAhqe GeV 2 
0.27±0M exp ±0.0A HQ E GeV 2 
0.20 ± 0.02 exp ± 0.02 H qe GeV 3 
-0.09 ± 0.04 exp ± 0.07 HQE GeV 3 . 



2.10 



> 

CD 


2.06 " 
















2.02 - 



0.10 
0.08 
0.06 
0.04 



"i 1 r 



(a) 



_l I L 



h 1 1 r 



_l I I L 




0.5 



1.0 



1.5 



1-2004 
8683A1 



> 17 
CD 

S 1.6 
=>< 

1.4 



( f ), 



0.5 



1.0 



(C) 



> 

CD 

CD 



0.18 



0.10 



1.5 



0.02 



(9), 



_l I I L 



0.5 



1.0 



-cut 



(GeV) 



"i 1 1 r 



> 

CD 



22 



20 



18 



"i 1 1 r 



-(d) 



0.00 
■S -0.01 
^-0.02 
-0.03 



"(h). 



1.5 



0.5 



_l I I L 



h 1 — "i r 



_l I I L 



1.0 



1.5 



Figure 1: Inclusive B — > X c ev decays: The measured hadronic-mass (a-d) and electron-energy (e-h) 
moments as a function of the cut-off energy, E cu t, compared with the result of the simultaneous 
fit (line), with the theoretical uncertainties [4] indicated as shaded bands. The solid data points 
mark the measurements included in the fit. The vertical bars indicate the experimental errors; in 
some cases they are comparable in size to the data points. Moments for different E cu t are highly 
correlated [8]. 



The fit results are fully compatible with independent estimates of /j>q = (0.35 ± 0.07) GeV 2 , 
based on the B* — B mass splitting [4], and of p\ s = (—0.15 ± 0.10) GeV 3 , from heavy-quark 
sum rules [9] . The lepton-energy and hadronic-mass moments have a different sensitivity to the fit 
parameters, but the results for the separate fits are fully compatible with each other and with the 
global fit to all moments. 

Moment measurements and the extraction of \V c b\ based on HQEs [10] were first performed by 
the CLEO collaboration [11]. More recently, global fits to a variety of moments were presented [12, 
13, 14] using HQEs in different mass schemes. The results presented here are compatible with these 
earlier fits, but have smaller uncertainties. 
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3 Exclusive B° 



D* l + v Decays 



The differential decay rate for the exclusive decay B — > D*~(. + v depends on three helicity am- 
plitudes which are commonly expressed in term of three form factors. For perfect heavy quark 
symmetry the three form factor have the same dependence on w and can be related in terms of the 
Isgur-Wise function. Here uu refers to the product of the four- velocities of the B and D* , which 
corresponds to the relativistic boost of the D* in the B rest frame, w = (Mg-t-M^,, — q 2 ) /2MbMd* . 

Exclusive B° — > D*~( + v decays are selected using primarily two kinematic variables, the mass 
difference Am = M^o _ — M—o and cos9b~d*£, where 9b-d*i represents the angle between the 
momenta of the B and the D*l pair. These distributions are also used to estimate the various 
background contributions from non-BB events, combinatoric and fake backgrounds to the signal 
D* and the leptons, and from decays involving higher mass charm mesons states. The backgrounds 
contribute less than 24% to the selected sample. 

For a sample of about 20,000 selected B° — > D*~e + v decays (restricted to D° — > K + ir~ decays) 
a binned maximum likelihood to the observed four-dimensional decay distribution was performed, 
see Fig. 2. This fit assumes a linear dependence of the form factors on w with a slope p 2 and form 
factor ratios, R\ and R2, that are independent of it;. 




1.0 1.2 1.4 -1.0 -0.5 0.0 0.5 1.0 -1.0 0.0 1.0 180 360 

w cosej cose v x 

8700A01 

Figure 2: Exclusive D* + e~v decays: Observed one-dimensional decay rates for the four variables a) 
w, b) cos 0£, and c) cos 0y, the helicity angles for leptons and hadrons, and d)x, the angle between 
the decay planes [15]. 



The results of the fit [15] 

^(^ = 1) = 1.328 ± 0.0m stat ± 0.025 syst 
R 2 (w = l) = 0.920 ± 0.048 stat ± OMSsyst 
p 2 = 0.759 ± 0.043 stai ± 0.032 syst 

represent a significant improvement over previous measurements [16]. 

The statistical errors still dominate; the largest contribution to the systematic errors is due 
to the background subtraction. Though the data do not allow for an independent determination 
of the u>-dependence of R\ and R2, they are compatible with predictions for this dependence by 
Caprini et al. [17] and others. 

To extract \V c h\ from exclusive B — > D*~i + v decays we integrate over the three decay angles 
and perform a least-squares fit to the observed w distribution. For this analysis we use 71,000 
decays, including electrons and muons and multiple decay modes of the D° meson. We rely on 
Heavy Quark Effective Theory (HQET) to relate the differential decays rate dr / dw to the product 
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Figure 3: Exclusive B° — > D*~l + v decays: Results of the fit to the decay rate: Top: the measured 
u; distribution compared to the fit result [18] (histogram). Below the signal events the various 
background contributions are indicated, A: D*X£v, B: uncorrelated D* — £ combinations, C: fake 
D decays, D: fake leptons, E: non-BB events , F: correlated D* — i combination. Bottom: the 
background- and efficiency-corrected data compared to the form factor parameterization with fitted 
parameters (solid line). 



\Vcb\ ■ A\ at w = 1. We have adopted the form factor parameterization [17] with 

Ri(w) ps i?i(l)-0.12(u>-l) + 0.05(w-l) 2 
R 2 (w) ~ R 2 (l) +0.11(u> - 1) - 0.06(w - l) 2 

and derive the w dependence in terms of a single unknown parameter, ■> 

AiH/ii(l)«l - 8p 2 Al z + (53p% -15)z 2 -{231p 2 Al -91)z 3 , 

where z = (VwTT - V2)/(Vw~+l + y/2) [17]. 

The data and the result of the fit are shown in Fig. 3. By extrapolation to w = 1 we extract [18] 

A!(l)\V cb \ = (35.5±0.3 siat ±1.6^)-KT 3 
p% = 1.29±0.03 stai ±0.27 syst . 

Using the lattice calculations for Ai(w) at zero-recoil of Ai(w = 1) = 0.919 ± 0.035 (including a 
0.7% QED correction) [19] we obtain 

\V cb \ = (38.7 ± 0.3 stat ± 1.7 syst ±H A1 ) ■ W- 3 . 

By integrating over the corrected w distribution we obtain for the branching fraction 

B D *-e+v = (4-90 ± Qm stat ± 0.36 syst )%. 

The dominant systematic error is due to the uncertainty in R\ and R 2 , for which we have used the 
earlier measurements by the CLEO Collaboration [16], i?i(l) = 1.18±0.32 and R 2 (l) = 0.71±0.21. 
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These values will be replaced in the future by the more precise BABAR results, taking into account 
also the w dependence of R\ and R2 including the non-linear terms. 

The measured CKM parameter \ V c b\ and the exclusive branching fraction B D *- e + u are consistent 
with earlier measurements based on the same parameterization [20] , except for those from the CLEO 
experiment [21]. 

4 Conclusions and Outlook 

With a significantly larger data sample and recent improvements in the theoretical calculation 
the BABAR Collaboration has succeeded in reducing the statistical and systematic error in the 
determination of The extraction of \V c b\ from exclusive decays has still sizable uncertainties. 
The current errors are dominated by the form factor uncertainties and the prediction of the decay 
rate at zero recoil. We expect further improvement, both experimentally and theoretical, in the 
near future. 
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